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ABSTRACT 



The perturbation method proposed by Professor Hui is 
described. The method gives exact solutions for the perturbed 
flow over both sides of a flat plate which is oscillating 
with small amplitude and frequency at large angles of attack 
in steady supersonic/hypersonic inviscid flow provided that 
the shock remains attached. Using the strip theory concepts 
these solutions are extended to study the dynamic stability 
in pitch of a flat, periodically oscillating wing or arbi- 
trary planform shape, at large angles of attack. Finally, 
Hui's perturbation method is extended to include the effects 
of upstream disturbances on a stationary wedge. 
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I. 



INTRODUCTION 



This thesis deals with oscillating bodies in inviscid, 
steady or unsteady, supersonic/hypersonic flow. Its sub- 
ject matter may be divided into three main topics. 

The first topic is covered in Sections II and III. A. 1,2 
and constitutes the background material for the other two. 

In Section II the Eulerian governing equations and boundary 
conditions for unsteady flow are formulated and the two- 
dimensional shock and expansion steady flow results are given. 
The basic elements of the linearized potential flow theory 
are also included in this section. In Sections III .A. 1,2 a 
perturbation method proposed by Professor Hui is presented 
[Ref. 6,7,8].^ The method uses as a basis the assumption 
that the unsteady flow over an oscillating flat plate, with 
attached shock waves at an arbitrary angle of attack, is a 
small perturbation from the steady reference flow and, for 
small amplitudes of periodic oscillations, it gives closed 
form solutions for the flowfield quantities in the disturbed 
flow regions. 

The second topic is covered in Sections III.A.3,B,C 
[Ref. 5], In Section III. A. 3 the closed form solutions over 
the upper and lower sides of the oscillating plate are combined 



Professor of Applied Mathematics, University of Waterloo, 
Ontario, Canada. NPS, Department of Aeronautics, Visiting 
Professor in the period January-August 1980. 
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and the strip theory approximation is employed to calculate 
the necessary quantities in the disturbed flow regions over 
a flat wing of arbitrary planform shape. The in-pitch stability 
derivatives for the three-dimensional wing are then obtained 
in closed form and the results are compared with other exist- 
ing theories. In Section III.C a comparison of these results 
with linearized potential flow theory results is included 
while, in Section III.B, the fundamentals of the linearized 
theory, as applied to three-dimensional wings, are presented. 

Finally, the third topic is covered in Section IV. Based 
on the same perturbation method this topic introduces upstream 
unsteadiness in the flow and its effects on stationary bodies 
are obtained. The upstream unsteadiness is of a fairly general 
periodic form to give the solution in the flowfield generated 
by a formation of bodies provided that the body originating 
shocks are not crossing and the expansion fan regions are 
not overlapping. An extension to the case of oscillating 
bodies is readily possible. 
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II. UNSTEADY INVISCID FLOW THEORY 



To describe the fluid motion two methods are available: 
the Lagrangian method and the Eulerian method. In both 
methods the fluid is regarded as a continuum, i.e., its 
matter is assumed to be continuously distributed. 

In the first method the fluid is assumed to be divided 
into infinitesimally small regions called fluid elements or 
fluid particles. The so-called particle point of view is 
then adopted and description of the fate of each individual 
fluid particle is sought . To determine the unknowns asso- 
ciated with each fluid particle, e.g., its position coor- 
dinates, density etc., a system of equations is set up by 
applying to each fluid particle natural laws such as Newton's 
second law of motion and conservation of mass and energy. 

These equations are known as the Lagrangian equations of 
fluid motion. 

Although the Lagrangian description appears to be a 
natural way to approach the problem of fluid motion, the 
Eulerian description is preferred in general since it gives 
more insight into the problem, it is much simpler and in 
most cases one is not interested in the fate of each individual 
fluid element but rather in the properties of the fluid at a 
certain point of the flow field at a certain time. 

In the Eulerian method attention is focused on the vari- 
ous points of the space filled by the flowing fluid and a 
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description of what is happening at each of these points, 
in terms of quantities of interest such as pressure, density 
and velocity is sought. The flow quantities of interest 
are assumed to be functions of space and time, i.e., to be 
scalar or vector fields. Thus in the so-called field point 
of view adopted here the fluid flow is characterized by the 
fields of velocity, pressure, density and so on and a fluid 
element occupying a certain point at a certain time assumes 
for its properties the values that are appropriate to that 
point at that time. To solve for these fields a system of 
equations is again set up by using, as before, natural laws 
such as Newton's second law of motion and conservation of 
mass and energy. These equations are known as the Eulerian 
equations of fluid motion. 

Throughout this thesis the Eulerian approach is used. 

In the next subsection the major steps in deriving the 
Eulerian equations of fluid motion are indicated and the 
equations are presented in the form in which they will be 
used later on. 

A. UNSTEADY EULER EQUATIONS 

The Eulerian equations of fluid motion, called hereafter 
simply equations of motion, may be set up either in differen- 
tial form or in integral form. Furthermore they may be 
developed either from the point of view of a certain fluid 
region that contains the same fluid elements for all times 
(control mass approach) or from the point of view of a fixed 
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volume in space through which different fluid elements flow 
through (control volume approach) . 

In what follows, the equations of motion are set up in 
the differential form from the point of view of an infinitesi- 
mal fluid region. 

In this derivation we will then naturally be involved with 
the calculation of the time rate of change of flow quantities 
as we follow the fluid element around, the so-called material 
derivatives of quantities. 

A physical interpretation of the material derivative’ and 
its components is briefly included. 

1 . Material Derivative 

Consider a fixed coordinate system and a fluid ele- 
ment situated at point r at some time t.^ Let Q(r,t) denote 
some fluid property Q of interest (density, velocity, etc.) 
associated with the point r at time t. The fluid element 
situated there (see Fig. 1.1) will assume for its corresponding 
property Q the value Q(r,t) . In a short time interval At 
the element moves through a distance As = VAt where V is 
its velocity at r and t. The element will then assume for 
Q the value appropriate to its new position r+VAt at time 
t+At. If we denote this value as Q (r+VAt, t + At) , the 
change of Q in the time interval At is 

AQ = Q (r + VAt , t + At) -Q(r,t) 



1 



Barred quantities denote vector quantities. 
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and the rate of change of Q following the element around, 
usually denoted by ~ is 



DQ = lim Q (r+VAt , t+At) - Q (r , t) 



Dt 



At-*-0 



At 



Expanding Q (r+VAt , t+At) in a Taylor series we get 



Q (r+VAt, t+At) = Q(r,t) + (f%_ At + (^-S)_ At 2 + 

C r , t 3t r , t 



+ (|§-)_ VAt + (2-S) (VAt) 2 + ... 
dS r , t as^ r , t 



where s denotes distance in the direction of the velocity 
V at point r and time t. 

Using the Taylor series expansion the rate of change of 
Q, becomes 



DQ dQ 9Q 

Dt 9t 3s 



( 1 . 1 ) 



where the derivatives are evaluated at r and t. 

The total rate of change of any property Q is thus com- 
posed of two parts . 

To see the physical interpretation of each of these terms 
consider a flow field which is at any instant spacewise uni- 
form but varies from instant to instant and a second flow 
field which is steady but not uniform spacewise. Consider 
also a fluid element in these fields which in a small time 
interval St is moving from position r in the flow field to 
position r+VSt. 
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V 




Fig. 1.1. Illustration of local, convective and 
material derivatives . 

The change of a property Q of the fluid element moving 

90 

in the first field described above is ^St and the rate at 

0 U 

which the property Q is changing locally at the point r is 

90 

This local rate of change is known as the local deriva- 
tive and is the first component of the total rate of change 
in equation (1.1). 

The change of a property Q of the fluid element moving 

90 

in the second flow field described above is (-r-=-)VSt. This 

dS 

change which is called convective change is necessary since 
the element has to have an appropriate value for its property 
Q at its new position r+V5t. The rate of change of the 
property is (^)V and is known as the convective derivative. 

d S 

It is the second component of the total rate of change in 
equation (1.1). 
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The sum of the local derivative and the convective 



derivative as given by equation (1.1) is known as the total 
or substantial or material derivative. The last term is 
probably more descriptive since the derivative is constructed 
following a certain material element around. This term will 
be used hereafter. 

Recalling that ^ represents the derivative of Q with 
respect to distance in the direction of the velocity V we can 
write equation (1.1) in the following form 

B¥ = It + <e v -9radQ)V 



where e^ is the unit vector in the direction of V or simply 



DQ 

Dt 



3Q 

at 



+ V*gradQ 



( 1 . 2 ) 



where Q can be a scalar or a vector quantity. 

If Q is vector quantity, i.e., Q = A the convective 
derivative V*gradA can be expanded using the formula 

V*gradA = j[grad(V*A) - V * curlA - A x curl V 

- curl(V x A) +V(divA) -A*(divV)] (1.3) 



2 . Equations of Motion 

We will now set up the basic equations that govern 
the unsteady motion of an inviscid, compressible fluid. We 
initially regard as unknowns the velocity field V(r,t), the 
pressure field p(r,t) and the density field p(r,t). We 
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want to establish relationships between these fields by 
applying to a certain fluid element the basic laws of nature; 
Newton's second law of motion, law of conservation of mass, 
law of conservation of energy. 



element situated at position r at time t (Fig. 1.2) . If 
V and p are the velocity and density of the element at r and 
t and if 6v denotes the volume of the element, then the mass 
and the momentum of the element are p6v and p6vV respectively. 
Let us also denote by F the total force acting on the element 
at time t . 



to any mechanical system states that "at any instant, the 

rate of change of momentum of a system is equal to the resultant 

of all forces that are acting on the system at that instant". 

By applying the above law to the fluid element considered 
and by noting that the rate of change of momentum of the ele- 
ment is simply the material derivative of the momentum we 
get 



a . 



Momentum Equation 



Let us consider an infinitesimally small fluid 



Newton's second law of motion which is applicable 



^r(p<SvV) 



F 



Since the mass of the element p6v remains constant the 



above equation becomes 




F 



(1.3.1) 
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-pndS 




d S = n d S 



Figure 1.2. Fluid Element 

The total force F is the resultant of the so- 
called surface forces and body forces. 

The body forces are forces that act throughout the 
body of the fluid such as the gravity force. These forces 
will be assumed to be small and will be neglected in this 
thesis . 

The surface forces on the other hand are internal 
forces in the nature of actions and reactions across the sur- 
fact that separates the fluid element from its neighboring 
fluid elements. For a frictionless or inviscid fluid the 
surface forces are simply pressure forces that act normal to 
the surface of the fluid element. Their resultant is 

- pndS 
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where p is the pressure of the element, 6S its total surface 
and n the outward unit normal to the surface. 

According to the integral definition of the gradi- 
ent of a scalar function the resultant of the pressure forces 
which is the total force F can be written 

F = - CfnS pndS = -6v gradp 

Equation (1.3.1) then becomes 

DV . ... . , 

Pi? - -gradp (1.4) 

This equation of motion is one of the fundamental equations 
of fluid dynamics and is called Euler's Equation. It repre- 
sents a system of three scalar equations for the five un- 
knowns — the pressure, the density and the three scalar com- 
ponents of velocity. 

b. Mass Equation 

Consider a fluid element situated at point r at 
time t with volume Sv (Fig. 1.3) . The mass of the element 
is pfiv. The law of conservation of mass states that "the 
mass of any fluid element remains constant as it moves about" 
even though, in general, its shape, volume and density may 
change . 

Applying the law of conservation of mass to the 
element considered is thus equivalent to setting the rate 
of change of mass equal to zero or setting the material 
derivative of mass equal to zero. This gives 
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0 



5t (pSv) = 



or 



-t Dp 
5V Dt + 



D 

Dt 



6 v = 



(1.5) 



The material derivative of the volume of the element may be 
expressed in terms of the velocity field. Let the surface 
of the element at time t be <5S and let it in a time inter- 
val <5t grow and become 6S^ as shown in (Fig. 1.3) . The 
change in volume of the element is equal to the volume swept 
by the surface of the element during the time 6t. If n is 
the outward normal to the original surface the net volume swept 
outward by 6S in time (St is given by 




Figure 1.3. Change in volume of a fluid element 
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V6 tnds 



where V is the velocity. 

The material derivative of the volume of the 
element is then if we also employ the definition of the 
divergence of a vector 



Dt 



Sv = 



V*ndS = 6v div V 



Equation (1.5) then becomes 



Dp 

Dt 



+ p div V = 0 



( 1 . 6 ) 



This equation is known as the equation of conservation of 
mass, or simply, the equation of mass, or the equation of 
continuity. It is a relation between the velocity and den- 
sity fields only. Since it does not involve any dynamical 
quantities (such as pressures or forces) it is a kinematical 
relation . 

c. Energy Equation 

The law of conservation of energy expresses the 
balance of energy exchanges that take place between a system 
and its surroundings. A fluid in motion may be regarded 
as a thermodynamic system characterized by the usual thermo- 
dynamic variables such as entropy, internal energy, etc. 

We will assume that the fluid is non-heat conducting and 
also that for a fluid element the only possible energy ex- 
change process is work done by the surface forces and body 
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forces. The law of conservation of energy applied to a 
fluid element may be expressed as follows. 

The rate of increase of energy E of a fluid 
element = the rate of work done by the surface forces and 
the rate of work W ^ done by the body forces. Symbolically 



DE 

Dt 



w x + 



W, 



( 1 . 7 ) 



We will assume, as before, that the body forces can be 
neglected and that the fluid is inviscid. The only possible 
energy exchange is then work done by the pressure forces . 

The rate of this work is 

W. = - <£jpndS*V » - <0 pV-ndS 

1 6 S 6 S 

where 6S as before, is the surface of the element, n is the 
outward unit normal and p, V are the pressure and velocity 
of the fluid element. 

According to the definition of the divergence of 
a vector we can write 

W, = pV • ndS = -6v div(pV) 

1 6S 

On the other hand the energy E of the fluid element 
is the sum of its kinetic energy and internal energy. We 
specify the internal energy of the fluid by the scalar field 
e(r,t) which denotes the internal energy per unit mass at 
point r and time t. Then since the kinetic energy per unit 
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mass is -j- the total energy of the fluid element with mass 
pSv is 

V 2 

E = p<5v(e + -—) 

Equation (1.7) can then be expressed as 

2 

P g^(e +^ 2 ") = - div(pV) (1.8) 

Equation (1.8) is referred to as the equation of conserva- 
tion of energy or simply the energy equation. An alternative 
form of this equation is 

P = - p div V (1.9) 

which can be found by subtracting from (1.8) the so-called 
equation of mechanical energy 

2 

P §t ( V } = " V'gradp 

The equation of mechanical energy is formed by multiplying 
both sides of Euler's equation (1.4) by V. 

The energy equation (1.8) has introduced the 
internal energy of the fluid as an additional unknown in the 
formulation of the governing equations. The list of unknowns 
includes the three scalar components of velocity, and also 
the pressure, density and internal energy of the fluid, while 
there are five equations available. 
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At this point we assume we are dealing with a 
perfect gas and introduce the equation of state for a per- 
fect gas 

p = p R T (1.10) 

as a sixth relation between the unknowns. 

We specify the temperature T of the fluid as a 
scalar field T(r,t) and since we are dealing with a perfect 
gas express the internal energy e of the fluid by the relation 

e = C v T (1.11) 



where C v is the specific heat at constant volume of the gas. 

From equations (1.9), (1.10) and (1.11) we can 

find that the energy equation for a perfect gas may be 
written in the following form which will be used hereafter 



D 

Dt 



Y 



= 0 



( 1 . 12 ) 



Summarizing we state that the basic equations that govern 
the unsteady motion of a non-heat conducting, inviscid, per- 
fect gas with constant specific heats are equations (1.4), 
(1.6) and (1.12). These equations are rewritten below for 
easy future reference. 



(continuity) 


+ V- (pV) = 0 


(la) 


(momentum) 


U + V-gradV + = 0 

3t 3 p 


(lb) 


(energy) 


3 (_P_) + v*v(-£-) = o 

3t \ y y 

P p 


tic) 

0 
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3 . Boundary Conditions 



Physical conditions that should be satisfied on 
given boundaries of the fluid are known as boundary condi- 
tions. There are several types of boundaries and consequently 
there are various possibilities for the boundary conditions . 

We will consider two types of boundaries which are 
of main importance in this thesis, Cl) "the solid-fluid 
boundary" where the fluid is bounded by a solid surface and 
(2) "the fluid-shock-fluid boundary" where two regions of the 
same fluid in different states of motion are separated by a 
flow discontinuity. The possibility of an infinitely weak 
discontinuity will not be excluded. 

The nature and number of the boundary conditions depend 
also on the form of the differential equations that govern 
the motion of the fluid. In this sense there are differences 
between the boundary conditions for a viscous fluid or an 
inviscid fluid. In the following the conditions for an invis- 
cid fluid are considered. 

a. Conditions at a Solid-fluid Boundary 

We assume that the fluid is bounded by an imper- 
meable solid wall and require that no fluid should cross the 
solid surface. Since the surface itself may be in motion we 
denote by V the velocity of the fluid and by V the velocity 

b 

of the surface. The relative velocity between the fluid and 
the surface is V - V . Let the equation of the surface be 
given by S(r,t) =0. A unit normal to the surface is then 
given by 
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— _ gradS 

± {gradS | 

and the component of the relative velocity normal to the 
surface is given by 

( V-V S ) -K = rIi i aF1 -[V-gradS -V s -gradS] 

Now assume that an observer moves with the surface particles 
that compose the solid surface. The observer cannot observe 
any change in the function S(r,t) considered as a scalar 
field. This means that the total rate of change of S(r,t), 
following a particle of the surface around, is zero, i.e., 

+ V *gradS = 0 

ot S 

The component of the relative velocity normal to the sur- 
face then becomes 



(V-V ) -n 
s 



± | gradF | 3 1 



(||r + V -gradS) 



(1.13) 



The condition of impermeability of the solid surface is 

(V-V ) -n = 0 

s 



or 



~ = — + V gradS =0 At S(r,t) = 0 (1.14) 

If the solid-fluid boundary is formed by the surface of a 
stationary rigid solid the above equation reduces to 
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V -grads 



0 



At S (r) 



0 



(1.15) 



It is pointed out here that the condition for- 
mulated above states that at each point of the solid-fluid 
boundary the normal to the surface component of the relative 
velocity between the fluid and the solid must vanish. Thus 
for an inviscid fluid nothing can be said about the tangen- 
tial to the surface component of the relative velocity which 
may or may not be zero. In short the so-called no-slip con- 
dition does not apply to an inviscid fluid. 

b. Conditions at the Fluid-Shock-Fluid Boundary 

Consider one dimensional adiabatic constant-area 
flow of a perfect gas through a discontinuity (Figure 1.4a) . 
Assume that the flow quantities in regions 1 and 2 are con- 
stant throughout the regions. The equations of continuity, 
momentum and energy between cross sections 1 and 2 give 
[Ref. Is pp. 55,56] 



P 1 U 1 = P 2 U 2 



Pi + 



P 2 + P 2 u 2 



u. 



u„ 



Y-l P- 



Y ! 
Y-l P- 



These equations hold as long as sections 1 and 2 are 
chosen outside the discontinuity region. The discontinuity 
region may be assumed to be vanishingly thin and sections 1 
and 2 may be brought arbitrarily close together. In this 
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Figure 1.4a. Change of conditions across a normal 
shock in a constant area duct. 

Steady flow. 




Figure 1.4b. Change of conditions across an arbitrary 

ciscontinuity surface S(r,t). Unsteady flow 
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case the requirement of a constant area duct is dropped and 
the equations apply locally across the discontinuity. Con- 
sidering the, discontinuity as a boundary separating regions 
1 and 2 we may refer to equations above as boundary conditions 
at the discontinuity and write symbolically 



[pu] = 0 Cl. 16) 

At x = x 

2 ° 

[pu + p] = 0 (1.17) 




Y 

Y-l 



£ 



0 



(1.18) 



where the square brackets denote the change in the enclosed 
quantity across the discontinuity. 

A generalization of the simple one-dimensional 
flow problem considered leads to the well known Rankine- 
Hugoniot conditions in the form that will be used in this 
thesis. Consider two regions in space separated by a surface 
S(r,t). Assume that adiabatic flow of an inviscid, perfect 
gas is established from region 1 to region 2 (Figure 1.4b). 
Assume also that the surface S(r,t) represents a discontinuity. 
In general the equation of the discontinuity is not known 
a priori but will be found as part of the solution of the 
flow problem. Thus the boundary is a so-called free boundary. 

Let the flow quantities in region 1 be V^(r,t), 
p^Fjt), p 1 (r,t) and in region 2 V 2 (r,t), p 2 (r,t), p 2 (r,t). 

The boundary conditions (1.16), (1.17) and (1.18) should apply 
locally at any point of the boundary S(r,t) provided that the 

o 
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velocities u^, u 2 are replaced by the normal components of 

velocities V , V . For a boundary moving with velocity 
n l n 2 

V g the normal to the boundary components of the relative 
velocities between the fluid and the free boundary are given 
by formula (1.13) 

\ - l l grW ff + ^ -gradS) 

V n 2 = ± Ti^dsT ( H + ' ? 2 ,gradS) 

From equations (1.16), (1.17) and (1.18) using the same 

notation we get: 

At S = 0 



(con tinuity) 


r /3S 
tp ( 8t 


+ V* VS) ] 


= 0 


(1.19) 


(normal 

momentum) 


r ,3S 

[p( Tt 


+ v-vs ) 2 + 


p(VS) 2 ] = 0 


(1.20) 


(energy) 


[i(^. +v 
L 2 v 3t 


•VS) 2 +-^- 
Y-l 


|(VS) 2 ] = 0 


(1.21) 



where, the square brackets denote, as before, the change in 
the enclosed quantity across the discontinuity and the symbol 
"V" stands for "gradient". 

The conservation of tangential momentum is not 
expressed by any of the equations above. To find the tan- 
gential momentum equation we require that the velocity com- 
ponent tangent to' the discontinuity be continuous. The 
tangential velocity components are given by 
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V = (n x V 1 ) x n 
*1 1 

V = (n x v_) x n 

where n is the outward unit normal to the surface 

n = VS = {S ,S , S } in Cartesian coordinates. The condition 
x y z 

of conservation of tangential momentum then becomes if the 
square bracket notation is used: 

(tangential [u] _ [v] _ [w] q = n 

momentum) S S S 

x y z 

where [u] , [v] , [w] denote the change of x, y and z components 

of velocity across the discontinuity. 

The equation above imposes two scalar conditions 
at the discontinuity, as one should expect from physical 
considerations. An alternative form of this equation is the 
following : 



(tangential 

momentum) 



[V x vs ] 



0 At S = 0 



( 1 . 22 ) 



Equations (1.19), (1.2)), (1.21) and (1.22) constitute the 

complete set of the Rankine-Hugoniot conditions. They are 
symmetrical and therefore remain unchanged if the brackets 
are taken to denote the change upstream rather than downstream 
through the discontinuity. A definite sense of flow direction 
is provided by the second law of thermodynamics , which requires 
that the entropy shall not decrease across a discontinuity. 

The change of entropy is given by [Ref. 1: p. 60]: 
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S -2p- = Int&VfY-D^-Y/fY-Dj 

K P]_ P]_ 

and the requirement stated gives the following condition 
which should accompany the Rankine-Hugoniot conditions 



(2nd law of 
thermo . ) 



[-^-] >_ 0 at S = 0 



(1.23) 



It should be noted that so far we have avoided using the 
term shock instead of discontinuity, though shocks are the 
only possible physical discontinuities. This was done on 
purpose since we intend to use the Rankine-Hugoniot condi- 
tions across hypothetical expansion fronts (negative or expan- 
sion shocks) through an iterative procedure so that in the 
limit condition (1.23) should not be violated. 

4 . Shock-Expansion Flows 

Unsteady flows with shock waves or expansion waves 
are considered in this thesis. These unsteady flows will be 
solved by first finding the corresponding steady flow solution 
and then using it as a reference flow in calculating the un- 
steady perturbation flow. 

Since a number of exact steady flow conditions are 
already available, they will be utilized as reference flows 
in finding the governing equations and boundary conditions 
of the corresponding unsteady flows. They include the super- 
sonic uniform wedge flow and the Prandtl -Meyer expansion flow. 
The results for these steady flows are stated below for easy 
reference. 
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a. Wedge Flow 



Consider steady uniform supersonic flow past a 



symmetrical two-dimensional wedge with semi-vertex angle 0 
(Figure 1.5). The wedge is assumed stationary. Oblique 
shocks will be formed at angles B measured from the free 
stream direction. By conservation of momentum the tangential 
component of velocity is continuous across the shock so that 
V. = V. . Then V and V are related by the normal shock 

H fc 2 n l n 2 

relations. Since V = V. sin 8 and V = V_ sin<t> the normal 

n ! 1 n 2 2 

shock relations can be used directly with replaced by 

sin B and replaced by sincf). The resulting relations 
for the oblique shocks are 




2yM^ sin^B - (y~1) 



Y+l 



« P | i 




Figure 1.5. Steady supersonic flow past a stationary 
wedge 
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2 . 2 
sin <j) = 



(y-l)M 2 sin 2 B + 2 
2yM 2 sin 2 3 - (y-1) 



(y+1) M 2 sin 2 B 
2 + (y-l)M 2 sin 2 3 



_2 _ tan 3 

>, tan <j> 



p 2 

From the last two relations eliminating — and recalling 



that <f> = 3-0 we get 

tan 8 = cot 3 



2 2 

M^sin 3 - 1 



y+t M 2 



.2 . 2 , 



2 - (M^sin 3 - 1) 



(1.24) 



One way to solve this equation, i.e., to find 3 for given 

and 0 is to express it as a cubic equation in x = cot 3 and 
select the appropriate positive root corresponding physically 
to the weak shock wave [Ref. 2: pp. 452-453], The following 
equivalent equation was used in this thesis for numerical 
calculations 



(1 + ^j^M 2 ) x 3 - (M 2 -l)cot9 x 2 +.(1 M 2 )x +cot 8 = 0 

where x = tan 3 . 

For attached shock waves this equation gives three 
real roots for 3 and the middle one is the one corresponding 
to the weak shock wave. 
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A maximum deflection angle x for given (beyond 
which the shock becomes detached ) can be found from equation 
(1.24) by differentiating it with respect to 6 and equating 
to zero . 



b. Prandtl-Meyer Flow 

Consider steady, two dimensional, uniform super- 
sonic flow over a convex corner (Figure 1.6a) . A turn of 
the flow through a single oblique expansion wave is not 
possible since this would lead to a decrease in entropy. 

The flow expands isen tropically through an infinite number 
of centered straight Mach lines that form the so-called 
Prandtl-Meyer expansion fan. Thus upstream of the ray OB 

where 0=0 the flow is uniform with Mach number M and 

00 

downstream of the ray OC where 0=0^ the flow is also uni- 
form with Mach number M^. For angles 0 such that 9^ <_ 0 <_ 0^ 
the flow field has the same properties along any ray 
0 = constant. The polar coordinate system shown has been 
chosen so that F = a + ^ + P (M^) where P (M) is the Prandtl- 
Meyer function 



P (M) = j- tan -1 [X (M 2 -l) 1/2 ] - tan" 1 (M 2 -l) 1/2 

with X = (^j) and y the ratio of specific heats of the 
gas. The Mach number M^ is given by 



For the wedge at an angle of attack a the flow deflec- 
tion angles are 0+a and 0-a for the lower and upper sides 
respectively . 
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Figure 1.6a. Steady supersonic flow over a convex corner 




surfoce 

fluid region 



Figure 1.6b. Expansion over a corner with a > I (M) x 
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P (M, ) = P (M ) + a 

1 00 



Along any ray 0 in the expansion fan the Mach number pressure 
and density are given by the following relations. 



0 

P 

P 



r tan -1 [A (M 2 -l) 1/2 ] 



= y /cy-d 

“ L E (M) J 

p t ^,l/(Y-l) 

P « l E(M) J 



where: 

E(M) = 1 + r2 



The r- and 0 -velocity components in the expansion fan are 
given by 



V = c sin 0 
r 

Vg = Ac cos 0 

where : 



c 



U 

00 




M £ 



It can be seen that the Prandtl-Meyer function 
P (M) defined above has a maximum value for M -*■ 00 . This value 
is 
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P (M) 



max 



J(/(Y+H/fy-l)-l) 



For convex corner angles a > P (M) 



max 



- P(M ) the streamlines behind 
00 



the expansion fan behave as though the flow occurred over an 

expansion of P (M) and for an inviscid fluid a region of 

max 

stagnant fluid lies between the hypothetical position of the 
body as sensed by the flow and the actual position of the 
body (Figure 1.6b). Whenever this occurred in numerical 
calculations performed, the influence of the flow on the sur- 
face of the corner^ was assumed negligible. 

B. LINEARIZED POTENTIAL EQUATION 

The so-called linearized theory of supersonic flow builds 
up the flow produced by the motion of a body by superposi- 
tion of small disturbances such as those produced by a moving 
sound source. One can develop in this way relatively simple 
methods for the computation of velocity and pressure distri- 
butions in the field. 

In the case of vortex-free flow, the equations of motion 
can be reduced to equations analogous to the wave equation. 
The coordinate parallel to the direction of the main flow 
plays the role of the time coordinate. Hence the methods 
of finding solutions of the wave equation can be used. 



We mean that in the equivalent case of a flat plate 
at an angle of attack a the pressure on its upper surface 
was assumed zero. 
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The linearized theory however has serious limitations. 
First, it gives only a first approximation since all devia- 
tions from the uniform parallel flow are considered small 
and therefore additive. This is justified for very thin or 
slender bodies at small angles of attack only. Second, there 
are speed ranges in which the linearization of the equation 
of motion even for small disturbances is not justified. For 
the linearized theory to be valid the following two conditions 
must be met. 

(a) The perturbation velocities must be small in compari- 
son to both the main stream velocity and the velocity of sound. 
This condition excludes the case of very high velocities 
since if the mean stream velocity is several times larger 

than the sound velocity, disturbances which are small rela- 
tive to the mean stream velocity may be of the same order of 
magnitude as the sound velocity. This speed range is called 
the hypersonic range. 

(b) The perturbation velocities must be small in comparison 
to the difference of the main stream velocity and the sound 
velocity. This condition excludes the range near M = 1, the 
so-called transonic range. 

In spite of the limitations described above the lineari- 
zation of the equations of inviscid, compressible fluids 
proved to be of excellent use in developing approximate solu- 
tions in the supersonic range. 

There are three general methods used in the linearized 
theory of supersonic aerodynamics . 
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(a) The method of fundamental solutions or sources. This 
method is based on the superposition of fundamental solutions 
of the linearized hyperbolic equation for small perturba- 
tions of a uniform supersonic flow. In formulating wing 
problems this method uses sources and doublets located in 
the plane of the wing and the strength of these singularities 
is determined so that the boundary conditions applicable to 
the wing planform and shape are satisfied. 

(b) The methods of acoustic analogy and operational cal- 
culus. In these methods the solution of the hyperbolic equa- 
tion is expressed by means of Fourier and Laplace integrals 
respectively. The second method is better adapted to super- 
sonic flow problems since the Laplace integrals exclude the 
possibility of upstream travelling signals while in the case 
of Fourier integrals one has to impose additional conditions 
to secure that this possibility is excluded. 

(c) The method of conical flows. This method is based on 
conical flows, i.e., flows for which the velocity components 
at points lying on a straight line drawn from a point chosen 
as vertex are independent of the distance from the vertex. 

In this method the solution of a hyperbolic equation in three 
variables is reduced to the solution of Laplace's equation or 
wave equation in two variables and the existing methods of 
conformal transformations and the theory of functions of' com- 
plex variables can be employed. 

Methods of higher approximations, i.e., methods which 
lead from the simple case of the linearized solution toward 

n 



42 



the exact one in successive steps, thus extending the range 
of satisfactory approximation, are required whenever the 
perturbation velocities are not small compared to the main 
stream velocity. 

In this section the concept of irrotationality will be 
introduced and then the basic assumptions and steps followed 
in deriving the linearized potential equation and the appli- 
cable boundary conditions for the general flow problem of a 
body in supersonic flow will be described. In the next sec- 
tion the linearized wing problem will be described. 

1 . Irrotational Flows 

By potential flow we mean that the velocity V is 
derivable from a scalar velocity potential <p , i.e., v = grad(j). 
On the other hand the vorticity or rotation oT in a fluid is 
defined as w = curl V and the flow is called irrotational 
if to = 0 or equivalently if 

3v _ 3u _ _3w _ 3v _ 3u _ 3w _ Q 

3x 3y - 3y 3z 3z 3x 

where u, v, w are the x-, y- and z- components of V. 

Physically the irrotationality of the fluid means 
that the fluid particles have zero moment; of .momentum about 
their own center-of— gravity ' axes or simply that they remain 
parallel to themselves as they move around. 

The condition of irrotationality of the flow is a 
necessary and sufficient one for the assumption of a poten- 
tial flow since the mathematical identities curl grad<£ = 0 
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and div curl V = 0 show that the velocity V can then be 
put in the form V = grade)) . 

The irrotationality throughout a flow field can be 
proved by applying the theorems of Kelvin and Stokes. 

Stokes' theorem states that the area sum of the rotation over 
a given area is equal to the integral of the velocity around 
a curve bounding the area. Formally, 

J to dA = j curl V • dA = <£> V*dX 
A 

The line integral is called circulation and is denoted by 
r. Thus F = 0 V*dX = Jto*dA. 

As a consequence of Stokes' theorem, to = 0 if the 
circulation r vanishes for all paths wholly within a simply 
connected flow region. 

Kelvin's theorem on the other hand states that the 

circulation T about any contour always composed of the same 

fluid particles (i.e., a fluid line) is constant in an inviscid 

fluid with only conservative or irrotational body forces . 

For an inviscid fluid it states that ^ = -0^- and it reduces 
DF 

to = 0 when there is a simple relation connecting p and 
p. Physically it means that circulation T about any line 
contour remains constant in time as we move along with the 
fluid. 

It is a consequence of both theorems that initially 
irrotational flows originating in a reservoir under uniform 
stagnation conditions or from straight parallel streamlines 
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will remain irrotational throughout the flow field at all 
times if there are no shock waves. 

2 . Linearized Potential Flow Equation 

In view of the above reasoning we assume that the whole 
flow field is irrotational and set V = grad$. As a result 
the three unknown velocity components u, v, w are expressed 
in terms of the scalar field 4> by 

u = $ , v = $ , w = 4 

x y ' z 

By using relation (1.3) the momentum equation (1.4) 
may be written as 

_ 2 

3V v 

p (^-£ + grad — - V x curl V) = -gradp 



or since curl V = 0 and v = grad$ we may write 



- » 



By integration we get 



3 $ Vi 

3t 2 



+ / 2P 
J P 



dp _ 



= F (t) 



If we define $ = $-/ F(t)dt this equation becomes 



* + + r = 

t 2 ■> o 



By differentiating this relation with respect to time 
and by taking its gradient we get the following two relations. 
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if the additional relation = gE- for the velocity of 
sound is used. 






2 

gradp 



We now introduce continuity equation (1.6) which, since 



Introducing in this equation the last two relations from 
momentum equation we get 



Equation (1.25) is the exact non-linear differential equation 
to be satisfied by the velocity potential $ for an unsteady, 
inviscid, irrotational flow. Because of its strong non- 
linearity, solutions have been found in very few special 
cases. Thus the small disturbance concept is introduced which 
leads to linearization of the equation. 



slightly in direction and magnitude from the free stream 
velocity , taken along the x-axis, and we define a distur- 
bance velocity potential V obtained from the total velocity 
potential $ by separating out the contribution of the uniform 
flow. In this way we have 



_ 2 

div V = V $, may be written as 



jr ft + l gradp + ' ,2i 



0 




(1.25) 



We assume that the velocity vector V differs only 
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Ux + ¥ 



4> = 

The local velocity components are then given by 

u = U + u ' , v = v 1 , w = w ' 

00 

where u' , v' , w' are the perturbation or disturbance veloci- 
ties which can be found from 



u 1 



= ’F 



x' 



V' = 



¥ , 

y 



w ' 



= ¥ 



We regard all disturbance velocities small in compari- 
son with U^, a and U^-a and all pressure and density changes 
small in comparison with main stream pressure and density. 

We also assume that the small quantities change gradually 
in all directions and that time variations are not too rapid. 

We now return to equation (1.25) and assume that 
the linear terms in Laplace's operator are of the same order 
of magnitude. Substituting the velocity V in this equation 
and retaining first order terms only we get 



v 2 t - 

a at 



+ 2U + U 2 i^l - 

00 8x3t °° « 2 

d X 



(1.26) 



where the second term in parenthesis is found from 



a 



=r[2$ $ J = 

2 x xt 



-Jj[2<D„ + u’)f xt ] = 



a 



a 



=r[2U ¥ . ] 

2 00 xt 



and the third term in parenthesis is found from 
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_L[<i> 2 $ ] = 

2 x xx 
a 



4uu„ + u')\ x ] ^ -4tu5 xx i 

a a 



To complete the linearization of (1.26) we should dispose 
of the factor 

a 2 

For steady flow a may be eliminated by using the 
following relation which is a consequence of the first law 
of thermodynamics 



2 , y-1 tt 2 , , 

a + ± ~ 2 ~ V = constant 

For unsteady flow let a slender body move in the flow 

field and denote the velocity of sound far upstream in the 

undisturbed field by the constant a m . In the vicinity of the 
2 

body a will then be a variable depending on position and 
can be represented by a sum of terms of the form 

2 2 . . .2 

a - a + (Aa) + ... 

00 



As a first approximation we set a = a m and introduce this 
value in equation (1.26) getting 



1 r Z 2 V 



2 

d V 



,2 d 2 v. 



^ + 2U_ ^hr+tT ] = 0 



~~ 2 L 7~T 

a co 



dxdt ’ Wc ° g x 2' 



(1.27) 



u 



An alternative form of this equation with M = — is 



co a 



2M . 

— V hr 'F =0 

" 08 ' "xx -yy zz a ro xt 2 tt 



(l-M )H' + 7 + Y 



(1.28) 



00 

Equation (1.28) is the linearized unsteady potential equa- 
tion and is used as the basic equation in most aerodynamic 
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analyses. Because of the assumptions made and terms retained 
this equation is valid for unsteady, inviscid, irrotational 
flows that are purely subsonic or purely supersonic and is 
limited to small disturbances only. 

3 . Linearized Boundary Conditions 

To completely specify the mathematical problem that 
describes the flow, the following boundary conditions need 
generally be prescribed [Ref. 3: pp. 1.27-1.29]. 

a) Surface Boundary Condition: The wing surface is 

impenetrable to the medium. 

b) Edge Conditions: Enough viscosity remains in the 

inviscid fluid to determine the flow pattern near sharp 
edges . 

c) Wake Conditions: The free vorticity shed from the 

trailing edge must have a circulation which vanishes together 
with the bound circulation. It is furthermore assumed that 
the shed wake is a continuous sheet of discontinuity which 

is coplanar with the wing projection in the direction of 
flight. Edge effects and rolling up of the sheet are 
disregarded . 

d) Conditions at Infinity: A state of uniform flow must 

be prescribed at inifinity. In addition the Sommerfeld 
radiation condition requires waves to propagate away from 
sources of disturbance toward infinity. 

e) Other Conditions: As the most important additional 

condition, the requirement that proper account be taken of 
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zones of influence and action at supersonic flow velocities, 
is mentioned. 

The first of the above conditions, namely, the con- 
dition of impermeability, as applied to a wing or airfoil, 
is considered next. The procedure can be easily extended to 
cover other bodies of interest such as slender fuselages, etc. 

Consider a wing fixed relative to a Cartesian co- 
ordinate system so that it lies close to the xy-plane (Figure 
1.7). Assume that the wing is submerged in an infinite mass 
of fluid moving with velocity in the positive x-axis. 

Let the upper and lower surfaces of the wing be expressed by 
equations 



S u = z - z u (x,y,t) = 0 
= z - z^(x,y,t) = 0 

The condition of impermeability of these surfaces, 
by recalling equation (1.14) , requires that 

3z 3z 3z 

w = TT + u Tx + v ~5y for z = V (x ' y) in R a 

3z i 3z i 3z i 

w = IT + u TTiT + v ~5y~ for z * V (x ' y) in R a 

where u, v, w are the components of velocity V and R is 

O. 

the portion of the xy-plane covered by the projection of the 
planform. 

These are exact, non-linear equations. To linearize 
them assume, as before, that the disturbance velocities are 
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Figure 1.7. Linearized boundary conditions for a wing 



small compared to the free stream velocity U and also that 
3 z 3 z 

the slopes — — , etc., are very small compared to unity. 

Then retaining first order terms only we get 



w 



w 



3 z 



u 



3z 



3 t 

3z * 
3 t 



+ U 



u 



00 3 x 

3z„ 



+ U 



3 x 



for z = z , (x,y) in R 

u 1 a 

for z = z^, (x,y) in 



Since and z^ are small compared to the wing chord 
we may, as a further step, replace the actual wing with an 
infinitesimally thick surface of discontinuities in u, v, w 
and pressure p. With this mathematical plane surface located 
on the xy-plane, we may expand w in Maclaurin series about 
its values just above and below the xy-plane 
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,+ 



w(x,y ,z u , t) 



= w 



<x,y,0 + ,t) +z„ *£ i xgj!> .. 'U + 



U 



w(x,y f z A# t) = 



w(x,y,0~,t) + z 3w<*-y-° ft) + 



Using the same arguments as before the higher order 
product terms can be neglected and the impermeability or 
flow tangency conditions take the following linear forms 



w(x,y, 




t) 



w(x,y,0 , t ) 
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u 



3 z 



3 1 



+ U 



u 



00 3x 



3t 



+ U 



!!i 

00 3x 



(x,y) in R. 



( x t y ) in R_ 



(1.29a) 

(1.29b) 
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Ill . OSCILLATING WINGS OF GENERAL PLANFORM 
IN SUPERSONIC/HYPERSONIC FLOW 



Consider a uniform (spacewise and timewise) supersonic 
or hypersonic flow of an inviscid perfect gas with constant 
specific heats past a flat wing of a general planform shape 
at an angle of attack a. Assume that the wing is performing 
a small amplitude slow pitching oscillation. 

The problem considered is to find the unsteady flow 
quantities in disturbed regions over the wing and thus its 
stiffness and damping derivatives. 

The governing equations of motion are given by equations 
(la-lc) , restated below. 

§| + V- (pV) = 0 (2.1a) 

|| + V-VV + ^2- = 0 (2.1b) 

ot P 

f-r(-^r) + V * V (~~) =0 (2.1c) 

p Y p Y 

where p, p, V and y are the pressure, density, velocity 
and adiabatic exponent of the gas. 

The flow tangency condition to be satisfied at the sur- 
face of the body is given by equation (1.14), restated below. 

|| + V-VS = 0 At S = 0 (2.2) 

< 3 1 

where S(r,t) = 0 is the equation of the body surface. 



u 
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The fluid-shock-fluid boundary conditions to be satis 
fied across the shock are given by equations (1.19-1.22), 
restated below. 



[ p<H + 


<1 

• 

<1 

£ 


= 0 At G = 0 


(2.3a) 




V* VG) 2 


+ p (VG) 2 ] =0 


(2.3b) 


[i(i£ + 

L 2 ' 3 1 


V* VG) 2 


+ : FIp (VG)2] = 0 


(2.3c) 


[V x VG] 


= 0 




(2.3d) 



where G(r,t) = 0 is the equation of the unknown shock shape 
and the square brackets denote the change in the enclosed 
quantities across the shock. 

Equations (2. 1-2. 3) are nonlinear. The nonlinearity of 
the governing equations and boundary conditions along with 
the existence of a shock with an unknown shape, contribute 
to the complexity of the problem considered. 

For low supersonic Mach numbers and very low angles of 
attack shock waves can be replaced by Mach waves and the 
linearized supersonic potential flow theory can be employed. 
The problem can then be solved, at least for certain groups 
of planform shapes, by fairly general methods. The fundamen- 
tals of the linearized supersonic potential flow theory as 
applied to a three-dimensional oscillating wing are presented 
in Section III.B. 

For high angles of attack and/or Mach numbers the shock 
waves become strong and the linearized theory cannot be used. 
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To overcome the difficulties encountered in this case, namely 
the existence of the shock and the nonlinearity of the equa- 
tions, Professor Hui proposed the use of a perturbation method 
in which the unsteady flow field is regarded as a small 
perturbation to some reference steady flow. Thus the unsteady 
flow problem is to be solved by first finding the corresponding 
steady flow solution and then using it as a reference flow 
in calculating the unsteady perturbation flow. The solution 
of the three-dimensional wing by this method will be presented 
in Section III. A. 

Finally, in Section III.C results for the stability deriva- 
tives are presented. A comparison with linearized potential 
flow theory results is also included. 

A. PROFESSOR HUI'S THEORY 

In this section the problem of dynamic stability of a 
flat wing of a general planform shape at arbitrary angles of 
attack in steady supersonic/hypersonic flow is considered. 

The wing is assumed to be oscillating in pitch with small 
amplitude and frequency and the bow shock be attached to the 
body at all times . 

The problem is covered in [Ref. 5] and only the basic 
steps will be included here, in Section III .A. 3. Its solu- 
tion is based on the assumption of an inviscid perfect gas 
with constant specific heats and the perturbation method 
developed by Professor Hui is employed to calculate the 
resulting unsteady flows over the upper and lower surfaces 
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of a two-dimensional flat plate. Finally, the strip theory 
approximation is utilized to combine the effects of these 
flows for the case of a three-dimensional wing, provided that 
the bow shock is attached and therefore the flows are independent. 

The unsteady flows over the lower and upper sides of a 
two-dimensional flat plate are studied in [Ref. 6,7] and 
[Ref. 8] respectively. Nevertheless, we will indicate, in 
the following first two sections, III.A.l and III. A. 2, the 
way in which these flow problems are formulated and solved. 

We will also give the solutions for the complete set of flow 
quantities in the disturbed regions. These flow quantities 
will be used in Section IV where the effects of upstream 
unsteadiness in the flow are considered. 

1 . Two-Dimensional Oscillating Flat Plate — Compression 
Side 

Instead of a flat plate the equivalent flow problem 
of a two-dimensional wedge is considered. This problem is 
formulated and solved in [Ref. 6,7] with the ultimate goal 
of studying the stability of wedges/caret wings. In what 
follows in this and the next subsection, 

a) The major steps in the method of solution are 
indicated . 

b) A generalized approach that permits the formulation of 
the fluid-shock-fluid boundary conditions is adopted. This 
approach is described in Appendix A and the formulation of 
the boundary conditions for the two cases is given in Appen- 
dices B and C. 
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c) Some of the results and discussions contained in 
(Ref. 6,7,8], pertaining to the stability of wedges/caret 
wings and flat plates, are not included since they are not 
directly related to the subject matter of this thesis. 

d) On the other hand the solutions for the complete set 
of the flow quantities in the disturbed regions, which are 
not included in the above references, are given. Much of 
the mathematical detail in obtaining these solutions is 
omitted. These flow quantities will be used, as mentioned 
before, in Section IV where the effects of upstream unsteadi- 
ness in the flow are considered. 

e) The same symbols as those used in the references will, 
in general, be employed. Changes will be limited- to those 
necessary for clarification purposes or generalization of 
approach. 

a. Problem Formulation 

Consider a two-dimensional wedge of length l, at 
design condition (zero mean angle of attack) , in a supersonic/ 
hypersonic, uniform, steady flow of an inviscid perfect gas 
with constant specific heats (Figure 2.1a) . Assume that the 
wedge is performing a low amplitude and frequency harmonic 
oscillation in pitch with given circular frequency w, about 
an axis perpendicular to the plane of the paper, through the 
point C shown. Let a system of cartesian coordinates Oxy 
be attached to the wedge so that 0 is at its apex and axis 
Ox is along the mean position of the upper surface. The bow 
shock is assumed to be attached to the body and the flow 
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wedge surface in 
steady flow 

un steady fb w 



shock surface in 

steady reference flow 
unsteady flow 



Figure 2.1a. Oscillating wedge in uniform steady flow 




plat e 



Figure 2.1b. Oscillating flat plate— Compression side 

u 
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quantities on the upper surface of the wedge are to be 
found. For all practical purposes the solution to this 
problem will give the flow quantities on the lower side of a 
two-dimensional flat plate of length Z sec 9 at an angle of 
attack 9 (Figure 2.1b) . 

b. Method of Solution 

The unsteady flow over the upper surface of the 
wedge will be found by perturbing the steady shock flow con- 
sidered in Section II. A. 4. a. 

Denote by U , p^, p ffi the velocity, pressure and 
density in region A. Denote by u q , p Q , p Q the velocity, 
pressure and density of the steady reference flow in region 
B. Non-dimensional lengths and time are introduced, defined 
by 

— — u 

x = ~ , y = and t = -3 t (2.4) 

Z Z l 

Assume that, as a result of the oscillation of 
the wedge, the perturbed flow quantities in region B are 
given by 

u = u + eu + ... (2.5a) 

o 

v = ev + . . . ( 2 . 5b) 

A 

p = p Q + cp + ... (2.5c) 

p’ =P o + ep + ... (2. 5d) 
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/\ 

where e describes the deviation of the unsteady flow from 
the reference steady flow. 

Substitute quantities (2.5) in the governing 
equations of motion (2.1), non-dimensionalize the indepen- 
dent variables, by using relations (2.4) and get the following 
perturbation equations 



u^ + u 
t x 



P u x 
o o 



V t + V 
t X 



p u l y 
o o 2 



p t + Px ■ 



p. + p + — (u +v ) = 0 

t K x u x y 
o 1 



(2.6a) 
(2.6b) 
(2.6c) 
( 2 . 6d) 



where subscripts denote partial differentiation and a Q is the 
speed of sound in the reference steady flow. Assume that 
the perturbation quantities u, v, p and p have the form 



u o e lkt U (x,y) 


(2.7a) 


u Q e lkt V (x,y) 


(2.7b) 


P o yM Q e ikt P (x,y) 


(2.7c) 


P o M Q e ikt R(x,y) 


( 2.7d) 



where U, V, P and R are unknown quantities to be found, M q 
is the Mach number in the reference steady flow and k is 
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the so-called reduced frequency of oscillation defined by 
k = a)2/u Q . k is assumed to be small and the time indepen- 
dent quantities U, V, P and R are expressed as power series 
in (ik) of the form 



U 



V 



P 



R 



(0) 

U v ' + 

v (o) + 

p'°> + 
R (0) + 



(ik)U (1) + ... 
(ik)V (1) + ... 
(ik)P (1) + ... 
(ik)R (1) + . . . 



(2.8a) 



(2.8b) 



(2.8c) 



(2 . 8d) 



Expressions (2.8) are substituted in (2.7) and the resulting 
expressions in the perturbation equations (2.6). By equating 
the same order terms, in each of these equations, a sequence 
of systems of partial differential equations is formed. 

Only the systems of zeroth and first-order equations are of 
interest in stability analysis. The zeroth-order equations 
are 



.(0) 


_L 


p (o) 


X 


M 

o 


X 


.(0) 


1 


p(0) 


X 


’ M o 


y 


(0) 


II 

O 




X 


X 




f (0> 


+ + 


M R ( 


X 


y 


O X 



(2.9a) 



(2.9b) 



(2.9c) 



(2 . 9d) 
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The first-order equations are 




-U 



( 0 ) 



(2.10a) 



V (1) + jjp- P ^ 
x M o y 




(2.10b) 




R ( 0 ) _ p (0) 



(2.10c) 





( 2 . lOd) 



Next we consider the boundary conditions applicable 
to the problem. 

Along the surface, the condition to be satisfied 
is the flow tangency condition given by equation (2.2). 

The equation of the surface for a stationary wedge is given 
by S = y = 0 and for an oscillating wedge is given by 
S(x,y,t) = y +e (h cos 0 - x) =0 where e = ee . Equation 
(2.2) then, gives with V = {u,v} = u q {1 + eU,eV}, 

V(x,y) = 1 + (ik)(x-h cos 0) at y = s(x-hcos0) 

Expanding V(x,y) about y = 0 and neglecting higher-order 
terms we get the linearized condition 

V(x,0) = 1 + (ik)fc-h cos 0 ) at y = 0 

Use of equation (2.8b) gives 




At y 



0 



( 2 . 9e) 
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V (1) = x - h cos 9 At y = 0 (2.10e) 

Across the shock, the conditions to be satisfied, are the 

conditions given by (2.3). The equation of the shock in 

the steady reference flow is given by G c = -y + x tan <p = 0 . 

Let the equation of the shock in the unsteady flow be given 

by G = -y + x tan tp + eQ(x) = 0, where Q(x) is an unknown 

function to be determined as part of the solution and 
A i. Jc t 

e = ee . To find the boundary conditions across the shock 
we substitute expressions (2.7) in (2.5) and the resulting 
expressions in equations (2.3). The boundary conditions, 
after linearization, are given by 



V = AQ' + (ik)BQ AT y = x tan <p (2.11a) 

P = CQ' + (ik)DQ (2.11b) 

U = EQ' + (ik)FQ (2.11c) 

R = GQ' + (ik)JQ (2. lid) 



where Q' = dQ/dx and the values for the constants A through 
J, which depend on the reference steady flow, are given in 
Appendix B. The derivation of these relations is lengthy 
and tedious, even for the case considered here, where there 
are no upstream disturbances. In Appendix A equations (2.3) 
are put in an alternate form. This form permits a much easier 

solution of the equations and is repeatedly used throughout 

o 

this thesis. 
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Expressing the time independent quantities, U, 
V, P and R by their expressions as given by equations (2.8) 
we get the following boundary conditions for the zeroth and 
first-order systems respectively. 



At y = x tan <p 


v (°> 


= 


AQ' 


(0) 


( 2 . 9 f ) 




p (°) 


= 


CQ' 


(0) 


(2 . 9g) 




u (0 > 


= 


EQ' 


(0) 


(2 . 9h) 




R (0 » 


= 


GQ' 


(0) 


(2 . 9i) 


at y = x tan <f> 


v<i> 


= 


AQ' 


(1 > + BQ (0) 


( 2 . lOf ) 




pU> 


= 


CQ' 


< X > + i Q (0) 


(2.10g) 




ad) 


= 


EQ' 


(i) + p Q (0) 


( 2 . lOh) 




R (1 > 




GQ' 


(1) + HQ (0) 


( 2 . lOi) 



c. Complete Solution 

Two boundary value problems have been set up. 

The zeroth-order equations (2.9a-2.9d) and the 
zeroth-order boundary conditions (2.9e-2.9i) constitute the 
first boundary value problem. This problem, which will be 
solved first, corresponds to the • problem of steady flow past 
a wedge and its solution should give the flow* quantities 
behind the shock for a stationary wedge. This result is 
shown in Appendix B. 



64 



The first-order equations ( 2 . 10a-2 . lOd) and the 
first-order boundary conditions ( 2 . 10e-2 . lOi) constitute 
the second boundary value problem. 

In both problems the equations and boundary con- 
ditions are linear and therefore suggest for the unknowns 
u(i), P ^ and R^, (i = 0,1), solutions that are 

linear combinations of the non-dimensional spatial coordinates 
x and y. 

In view of the above reasoning, we assume for the 
first problem a solution of the form 



U 



( 0 ) _ 



( 0 ) 

‘1 



u ' x + u 1 "" y + u 



(0) 

‘2 



( 0 ) 

3 



(2.12a) 



V 



( 0 ) ( 0 ) , ( 0 ) . ( 0 ) 
= v£ x + Y + v^ 



( 0 ) _ ( 0 ) 



R 



Pl x + p 2 y + P 3 



( 0 ) _ ( 0 ) 



r i x + r 2 y + r 3 



( 0 ) _ ( 0 ) 



q l X + q 2 



( 0 ), 

2 ■ 

( 0 ). 

2 

( 0 ) 



( 0 ) 



( 0 ) 



(2.12b) 
(2.12c) 
( 2 . 12d) 
(2 . 12e) 



Substitution of these values into the zeroth-order equations 
and boundary conditions (2.9) gives the following solution 



U 



V 



(0) 


o 

^ CO 

3 

11 


= E/A 


(2.13a) 


(0) 


(0) 

= V 

3 


= 1 


(2.13b) 


(0) 


- / 0) 


= C/A 


(2.13c) 
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G/A 



x/A + q 



( 0 ) 

2 



x/A - h cos 8 



(2.13d) 
( 2 . 13e) 



where the coefficient was found from the additional con- 

dition that the bow shock is attached to the leading edge 
of the body. Recalling that the equations of the surface 
and the shock are given by S(x,y,t) = -y+e (x-hcos0) = 0 
and G(x,y,t) = -y+xtan<J>+eQ (x) = 0 we get the above result 
by letting x = 0 in the equation e(x-hcos0) = xtan4>+s(Q^ 

+ (ik)Q (1) ) . 



Similarly to solve the second boundary value 
problem we assume a solution of the form 



(1) 


= u/’x ♦ 


u/’y 


♦ u' 1 ’ 


(2.14a) 


(1) 


= v/'x ♦ 


v/’y 


+ v/» 


(2.14b) 


(1) 


= p‘ ll x + 


p/’y 


+ p/’ 


(2.14c) 


(1) 


= r/>x + 




+ r/> 


(2 . 14d) 


(1) 






X ♦ q/’ 


(2 . 14e) 



and substitute these values to the first-order equations and 
boundary conditions (2.10). The resulting solution is given 
below 



u^ = - [Cq| 1) /M + (E +D/M q ) /A + 2tan<j> ] 
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u 



u 



V 



V 



( 1 ) 

2 

( 1 ) 

3 

( 1 ) 

1 

( 1 ) 

2 

>[ 1! 

■l 11 

>/> 

.( 1 ) 

'1 

.( 1 ) 

2 

.( 1 ) 



( 1 ) 



[qj 1) (E+C/M Q ) + (E+F+D/M o )/A]/tan <f» +2 
h COS 0 [E (B-l) /A - F] 

1, v 3^ = “ h cos 9 
(Aq ; [ 1) +B/A - 1) /tan <p 
Cq^ 1) + D/A + 2 M q tan <p 
- 2 M 

o 

h cos 9 [C (B-l) /A -D] 

Cqj 1) + (D +C -G)/A + 2 M q tan <p 
t(G - C) q^ 1 ^ + (G + H -C - D) /A] / tan <f>- 2 M q 
h cos 0 [G (B-l) /A - H] 

[k 2 (A-B-C+E/M q ) /A -M tan <p[D/A+2 M tanf) ]/2(k 2 A 

+ CM tan<}>) 
o 



2 2 2 
k = M7(M -1) 
o' o 

2 . Two-Dimensional Oscillating Flat Plate — Expansion 
Side 

This problem is formulated and solved in [Ref. 8] . 
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a. Problem Formulation 



Consider a two-dimensional flat plate of length 
I in a supersonic/hypersonic uniform steady flow of an 
inviscid perfect gas with constant specific heats. Assume 
that the plate is performing a low amplitude and frequency 
harmonic oscillation about its leading edge. Let a system 
of Cartesian coordinates Oxy be attached to the body so that 
0 coincides with the leading edge of the plate and axis. Ox 
is along its mean position (Figure 2.2a). The bow shock is 
assumed to be attached to the leading edge and the flow 
quantities on the upper surface of the plate (in Region C) 
are required. 

b. Method of Solution 

The unsteady flow, over the upper surface of 
the oscillating plate, will be found by perturbing the 
steady Prandtl-Meyer flow, considered in Section II.A.4.b. 
Denote by U^, and the velocity, pressure and density 
in Region A. Denote by u^, p^, the velocity, pressure 
and density of the reference steady flow in region C. The 
solution procedure that follows is similar to the one pre- 
sented in Section III.A.l.b and most of the assumptions and 
results given there, apply to this section too, provided 

that u , p , p , M and a (the reference steady flow quanti- 
o *o o o o 

ties over the compression side of the flat plate) are replaced 
by u^, p 1# p , M 1 and (the reference steady flow quantities 
over the expansion side of the flat plate). 

o 
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Figure 2.2a. Oscillating flat plate — Expansion side 
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o 




Figure 2.2b. Oscillating flat plate--Polar coordinates 
for P rand tl -Meyer flow 



Non-dimensional time and lengths are introduced, 
defined by (2.4). As a result of the oscillation of the 
plate assume that the perturbed flow quantities in region C 
are given by (2.5). The resulting perturbation equations are 
given by (2.6). Assume that the perturbation quantities are 
given by (2.7) and let the time independent quantities U, V, 

P and R be expressed as power series in (ik) by (2.8) . The 
zeroth and first-order equations, derived as before, are 
given by (2.9a-2.9d) and ( 2 . 10a-2 . lOd) restated below. 
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- . _L p<°) 
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M l x 


-(0) 




= - p(°) 
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i — 1 
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(0) 




= R <°> 
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V <0 > + M, R (0) = 
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y lx 
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II 
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-(1) 
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P 

II 
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X 
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(1) 




R (1) = R (0) - P ( 


X 




X 


f U) 


+ 


V (1 > + M. R (1) = 


X 




y lx 



- M. R 



( 0 ) 



(2.16a) 
(2.16b) 
(2.16c) 
( 2 . 16d) 
(2.17a) 
(2.17b) 
(2.17c) 
( 2 . 17d) 



We now consider the boundary conditions applica- 
ble to the problem. 
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Let the equation of the oscillating surface of 
the plate be given by S(x,y,t) = y - ex = 0.^" The flow 

tangency condition (2.2), with V = u 1 (l+£U,eV} / gives after 

linearization 

V (0) =1 At y = 0 ( 2 . 16e) 

V (1) = x At y = 0 ( 2 . 17e) 

We next consider the boundary conditions across 
the surface that is separating regions B anc C. We assume 
that as a result of the small amplitude slow oscillations 
of the body the separating surface is slightly deformed and 
its equation is given by G = -y + x tan tj) + eQ(x) = 0 where 
Q (x) is an unknown function which may be expanded as 
Q(x) = + (ik)Q^ + ... . We call the flow expansion 

an expansion front or, simply, front and assume 
that upstream of it the Prandtl-Meyer flow is not disturbed 
while, along the front, the unsteady flow matches the steady 
Prandtl Meyer flow continuously. The assumptions made are 
completely analogous to the assumptions made in the case of 
a finite compression shock discontinuity and the Rankine- 
Hugoniot conditions (2.3) may be used to give the boundary 
conditions across the expansion front. We note that, since 



The difference in the form of this equation and the one 
considered in the previous section is due to our assumption 
that in this case the flat plate is oscillating about its 
L E and thus h = 0 . 
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we are dealing with an infinitely weak discontinuity, the 
flow across the front is isentropic. The procedure is given in 
Appendix C and the resulting boundary conditions, after lineari- 
zation and' use of equations (2.8), are given below. 



V 



V 



U 



(0) 


= 


A'Q’ (0) 


(0) 


= 


C ' Q ' (0) 


(0) 


= 


E'Q' ^ 


(0) 


= 


G'Q' (0) 


(1) 


= 


A'Q' (1) 


(1) 


= 


C ' Q ' (1) 


(1) 


= 


E'Q' ^ 


(1) 


__ 


G'Q' ^ 



At y = x tan <j> 



At y = x tan <J> 



’Q 

'Q 

'Q 



( 2 . 16f ) 
( 2 . 16g ) 
( 2 . 16h) 
( 2 . 16i) 
(2 . 1 7 f ) 
( 2 . 17g) 
(2 . 17h) 
( 2 . 17i) 



where the coefficients A' through J' depend on the reference 
steady flow and are given in Appendix c . 
c. Complete Solution 

Two boundary value problems have, again, been 
set up and will be solved successively. The zeroth-order 
equations and boundary conditions (2.16) constitute the 
first problem, which corresponds to a steady 
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Prandtl-Meyer 


flow problem 


. To solve 


it 


assume , 


as before. 


that the unknowns 


have 


the 


form given 


by 


(2.12) , 


substitute 


in (2.16) and 


get 


with 


K = 


m l / >/mJ - 1 








p(0) 


= 


p< 0) 


s 


K 






(2.18a) 


R (0 > 


= 




= 


< 






(2.18b) 


u<°> 


= 


(0) 

U 3 


= 


-</M 1 






(2.18c) 


V<°> 


= 


(0) 

V 3 


= 


1 






( 2 . 18d) 


Q (0) 


=5 


r — 1 

+ ^ 

>- 


4 

X 


+ g (0) 
q 






( 2 . 18e) 



where the coefficient = 0 since, in this case, the 

plate is oscillating about its LE . 

The second boundary value problem consists of 
the first-order equations and boundary conditions given by 
(2.17). In this case the unknowns are expressed by (2.14) 
and the resulting solution is 
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3 . Three-Dimensional Oscillating Wings of Arbitrary 
Planform Shape 

The problem is formulated and solved in [Ref. 5] 
and only the major steps in the method of solution are given 
here . 



a. Problem Formulation 

Consider an oscillating wing of arbitrary plan- 
form shape at an angle of attack in a steady, uniform, super- 
sonic/hypersonic flow (Figure 2.3). Assume that the oscilla- 
tions are periodic with small amplitude and frequency and 
that the bow shock is attached to the wing. We let the 
pressure, density, velocity and Mach number of the approach- 
ing flow be given by p^, p^, and M^. We also denote the 
total area and the root chord of the wing by S and Z respec- 
tively and the distance of the pivot position from the leading 

edge by x . We assume that the pitching motion of the wing 
c 

is described by 



6 (t) 



0 



iwt 

e 



where 9* and w are the amplitude and circular frequency of 
oscillation and t is the non-dimensional time. 

We define the reduced frequency of oscillation 
by 



k = u)J£/U 

' oo 
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b. Method of Solution 



The in-pitch stability derivatives of the wing 

are required. The pitching moment coefficient C^, the 

stiffness derivative -C and the damping in-pitch deriva- 

m Q 

tive -C are defined by 
m (3 

= 2M/p oo U 2 £S = (1/JtS) / / (x-x c )C p (x,y,t)dS 

- e (t) [ (-C ) + (ik) (-C ) ] (2.20) 

m Q 



where M is the moment about the pivot axis and C is the 

P 

pressure coefficient. 

The pressure coefficient is defined as usual by 



2(P ~Pj 

P U 2 
^00 00 



and using Idle two-dimensional flow assumption can also be written in the form 



x 

C P = (C P } + 0 (t) f A + ik( B J -C -£) ] 

= (C ) + 0 (t) [A +i^(Bx -Cx c ) ] (2.21) 

where A, B and C are dimensionless constants which are 
functions of geometry and steady flow quantities and (C ) Q 
is the mean pressure coefficient corresponding to 0(t) = 0. 

This coefficient has no contribution to the stability 
derivatives and will be neglected. 

Our goal is to calculate the coefficients A, B, 

C appearing in relation (2.21) and use the pressure coefficient 
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in relation (2.20) to obtain the stability derivatives. 

It is clear that, for a two-dimensional flat plate with 
attached shock wave, the coefficients will be the sum of 
independently obtained coefficients over the compression 
and expansion sides of the plate respectively. Thus the 
coefficients will have the form 

A - V B - B t - V C = C ! - C u <2 - 22 > 



where the subscripts £ and u denote the lower and upper 
surface respectively. 

The results from Section III.A.l and [Ref. 6,7] 
are used to obtain coefficients , B ? and of the form 

A Z = A q C/A, B £ = y Q (2G-I), = U o I (2.23a) 



where X and y are coefficients introduced to account for 
o o 



differences in notation given by 



- p u - 
(- 2 .) 2 

M M U ' ’ 

O 00 oo 



y = X U /u 
o o °° o 



and G, I are quantities defined in [Ref. 6,7]. 

Similarly the results from Section III. A. 2 and 
[Ref. 8] are used to obtain coefficients A^, B^ and of 
the form 



A u - -W<M 2 -1)°- 5 , C u - Vl /X l 

B u = u l M l( M l'2 , / (M l‘ 1)1 ' 5 ' 
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where X^, y^^ are correction coefficients to match present 
notation of the form 



2 P 1 U 1 2 

557<^> <l£ ■ 



= X.U /u. 
1 00 1 



c. Solution 

Substituting relations (2.23) in relation (2.22) 
and putting the resulting expressions in relation (2.21) 
we obtain the pressure coefficient for a two-dimensional 
flat plate. We shall employ the strip theory to solve the 
three-dimensional wing problem in hand, which means we will 
assume that the flow, at each point of the wing, is two- 
dimensional locally. This assumption permits the use of the 
two-dimensional flat plate pressure coefficient for the case 
of the three-dimensional flat wing. 

Using relations (2.21) and (2.20) we get 



C 

m. 

9 



where 



-C = A (I. - x /%) (2.24a) 

m 9 1 c 

[bi 2 +(c-b)i 4 ]-[(b+c)i 1 +(c -b)i 3 ]x c /ji+ C(x c /Ji) 2 

(2.24b) 



I, = k / (g 2 -f 2 )dn, I 9 = 2k/3 / (g 3 -f 3 )dn (2.25) 

1 0 ^ 0 

1 1 2 2 
I. = 2k. / f(g-f)dn, I, = k / f(g -f ) dn 

J 1 0 0 

k = Vo/s , u n = y/b 
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For delta wings with power law leading edges, i.e., wings 



delta wings may be obtained, which is practically independent 

of the power n, as shown in (Figure 2.4). Plots of C and 

m 9 

C vs a for several values of pivot axis position and powers 

m 9 

n are given in (Figure 2.5). In (Figures 2.6, 2.7) the 
stability derivatives vs the pivot position are plotted for 
several values of power n and angles of attack 10° and 20°. 
Comparisons of results obtained by the present theory with 
results obtained by other theories and related discussions 
are included in [Ref. 5] and will not be repeated here. 
Nevertheless in Section III.C we will compare the present 
results with potential flow theory results. 



with f(ri) = and g(n) = 1, equations (2.25) become 



I 



1 



(n+1) / (n+2) , 



I 



2 



(n+1) / (n+3) 



(2.25a) 



I 



3 



n/ (n+2) 



I 



4 



n/ (n+3) 




x /% = [B+ (2n+l)C]/2C(n+2) 



is 




min 



1 r B+ ( 2n+l) C ,2 
C 1 n+2 J 



•] + 4 (B+nC ) / (n+3) 




= 0 a stability boundary for power law 
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Figure 2.4. Stability boundary 
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Figure 2.5. Stability derivatives vs AOA 
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* , 

curves continued to line through 

*c/l= I. 0 

Figure 2.6. Stiffness derivative vs pivot position 
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Figure 2.7. Damping derivative vs pivot position 
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B. LINEARIZED THEORY 



It was pointed out in Section II that several methods 
are used in the linearized theory of supersonic aerodynamics. 
A simple application of one of these methods, the so-called 
method of fundamental solutions, will be described in this 
subsection. For this application a distribution of pulsating 
sources over a "simple" planform will be considered. 

Before dealing with this application, however, the notion 
of a disturbance propagation in supersonic flow will be re- 
viewed and the fundamental solution for a moving pulsating 
source will be introduced. 

1 . Propagation of Disturbances 

Supersonic flow is dominated by the fact that dis- 
turbances travel with finite velocities, namely, the speed of 
sound. In formulating the linearized potential equation 
(1.28) deviations of the speed of sound from its free stream 
value a m were neglected and also the perturbation velocities 
u' , v' , w' were taken very small compared to the free stream 
velocity U^. As a result any disturbance from a source 
located at point (x,y,z), in a coordinate system fixed to 
a body in the flow, can be felt only inside or on the surface 
of a right cirular cone whose axis points downstream from 
the source (Figure 2.8). An observer moving with the fluid 
sees a pulse emitted at t = 0, expanding on a spherical sur- 
face with instantaneous radius a t and center moving down- 

CO 

stream with velocity U^. The positions of these expanding, 
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Upstream zone of influence Downstream zone of influence 




Figure 2.8. Upstream and downstream zones of influence. 

Disturbance located at (£, n,?) at time t = 0 

moving spheres form an envelope which is a cone of semivertex 

angle 

-1 “ oo 1 " .-11 -11 

y = sm - = sin — = tan 

U t M /— = — 

00 oo / 2 

VM -1 
00 

where y is the Mach angle. 

This cone is known as the Mach cone or downstream 
zone of influence of the point (x,y,z) and its equation is 
given by 



(C-x) 2 - (M^-l) [ (n-y) 2 + U-z) 2 ] = 0 (2.26) 

On the other hand the point (x,y,z) can be influenced 
by sources whose locus is evidently a similar cone directed 
forward from (x,y,z) . This cone is known as the forecone 
or upstream zone of influence and has the same equation except 
that for this case x > E, . 
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Assume now that a steady supersonic flow has been 
established past a stationary three dimensional wing lying 
very close to the xy-plane. Regarding each point of the wing 
as a disturbance source we can see that the downstream zone 
of influence for the entire wing is bounded by the envelope 
of the Mach cones emanating from the leading edge. If the 
leading edge is straight (or the wing is two-dimensional) 
the envelope reduces to the so-called Mach wedges . 

To calculate the fluid motion at any point (x,y,z) 
we need to consider only the contribution from the disturbance 
sources that belong to the region of the xy-plane intercepted 
by the forecone from point (x,y,z). This area of influence 
forms a hyperbola and, with sources assumed to lie on the 
xy-plane, is found from the equation of the forecone (2.26) 
by setting 5 = 0. We thus get 

n 1/2 = y ± j (x-o 2 /(M^-i) -z 2 (2.27) 

A generalized supersonic planform with leading edge 
AA'C'C, trailing edge DD'F'F and streamwise tips AF and CD 
is considered next (Figure 2.9) . For each point (x,y,0 + ) the 
area of the sources that influence the point reduces to a 
region bounded by two straight lines upstream of the point. 
These lines found from the forecone equation (2.26) by setting 
z = 0 and 5=0 are given by 

£-x = ± vtr-i (y-n) £ < x 

oo — 
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\ 



Figure 2.9. Generalized supersonic planform 

and are shown as dotted lines making an angle y with the 
x-axis (Figure 2.9). 

The portions A'BC' and D'EF' of the leading and 
trailing edges are called supersonic since the velocity normal 
to these edges is greater than the speed of sound. Similarly 
the remaining portions of the leading and trailing edges are 
called subsonic since the normal component of velocity is 
less than the speed of sound. Along the supersonic portions 
of the leading and trailing edges there is no communication 
between the upper and lower surfaces of the wing and the flow 
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over the top or bottom of the planform can be calculated 
without reference to the shape of the opposite side . On 
the contrary the upper and lower surfaces of the wing are 
not independent along subsonic portions of the edges . 



to illustrate different sorts of upstream influence regions. 
The difficulty in solving the linearized flow problem, i.e., 
finding the fluid motion at these points, increases as we 
move from point 1 to point 3. 



There exists no universal method of approach in solving 



the linearized problem for different sorts of influence re- 
gions. Thus each planform shape calls for a different method 
of approach. 



supersonic leading and trailing edges, the so-called simple 
planform, that will be considered later in this section. 

2 . Fundamental Solution of a Moving Source 

For = 0 the linearized potential equation (1.28) 
reduces to 



This is the classical wave equation for the propagation of 
sound in a still medium and its fundamental solution is given 
by 



Points 1, 2 and 3 shown (Figure 2.9) are selected 



It is the simplest case of a planform with purely 



y + y + y 

xx yy zz 




T(x,y,z,t) = ^ F ( t — --) 

00 



(2.28) 
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where is the speed of sound of the undisturbed flow and 

• . / 2 2 2 
r is the radial distance from the origin, i.e., r = vx +y +z . 

A solution of the linearized potential equation can 
thus be readily found if we can transform it into the classi- 
cal wave equation. To achieve this the following transfor- 
mation, known as the Lorentz transformation, is employed. 
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T 



\/l -M 2 



x, n = y y, 



ia vl - M 2 

OQ 00 



•(t + 



M 



a ( 1 - M ) 

OO 00 



x) 



y z 



The velocity potential of a sound source fixed with 
respect to the xyz system is then found to be 



’F (x,y , z , t) 



F [- 



ya vi^-M 2 



A -M 2 



R 



-M x + R 

OO 00 

-(t )] 



a ( 1 - M ) 
00 00 



f~2 2 2 2 " 2 

where R = v x + (1-M ) (y +z ) and 1-M is known as the 

00 - 1 00 

Prandtl factor. 

To find the constants a and y introduced by the 
Lorentz transformation we require that the sound source should 
produce constant sound flux independent of the free stream 
Mach number and we get 



a 



a .* 1 



y 



A -M 2 

00 



Thus the velocity potential of the moving source becomes 
if in addition the relation A = - 7 ^- is used. 
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¥ (x,y , z , t) 



4ttR 



F(t - 



-M x + R 



a (1 - M) 
00 00 ' 



Expressing the above result in a coordinate system moving 
uniformly with the source located at (£,n,?) we obtain 



f(x,y,z,t) = -Ji5 p( t-r-) “ - 4 iR F(t - T l> (2 - 29) 

00 

where : 

R = / (x-7) 2 + (1 -M^) [ (y -n) 2 + (z - ~) 2 ] 

and 



T _ D _ -M(x-g) + R 

^ a °° a (1 -M^) 

00 00 

By comparison of expressions (2.28) and (2.29) it is seen 
that the solution for a moving source can be obtained from 
the solution for a stationary source by replacing the ordinary 
distance r by R in the amplitude and by D in the phase. The 

quantities R and D are called amplitude and phase radii 

respectively. 

A geometric interpretation is given in Figure 2.10. 

At time t a field point Q and a source 0 moving with velocity 
in the negative x-direction are considered. For supersonic 
flow there are two spherical waves passing through point Q 
at time t. These waves originated from the source at times 
t- and t-T^f at which times the source was located at 
positions and shown (Figure 2.10a). For subsonic flow 
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R s 0 




Figure 2.10a. Supersonic source 




Figure 2.10b. 



Subsonic source 

o 
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there is one spherical wave passing through point Q that 
originated at time t-T when the source was located at point 
P (Figure 2.10b). From the geometry of the figure we get 

(x-c-u T) 2 + (y-n) 2 + (z-O 2 = a 2 T 2 (2.30) 

Solving this equation for time T we get for subsonic flow 
one real positive solution and for supersonic flow two real 
positive solutions. 

Physically we are looking for the effect that a dis- 
turbance, originating at point (£, n,£) at some time t-T, 
will have at some later time t at a point (x,y,z). In this 
sense the potential is a retarded potential. For supersonic 
flow the disturbance is first felt at some point (x,y,z) 
after a certain time has elapsed. The point (x,y,z) 
penetrates the wave front of the disturbed region and be- 
cause. it is moving at a speed greater than that of the wave 
front it emerges from the disturbed region at some later time 
• For subsonic flow once the point (x,y,z) penetrates the 
wave front (after a certain time T has elapsed) it will remain 
in the disturbed region since its speed is less than that of 
the wave front. Finally the nonexistence of positive roots 
of equation (2.30) should be associated with an undisturbed 
region, i.e., with 'F = 0 . 

In view of above reasoning the source solution for 
supersonic flow takes the form 
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